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Abstract: In this paper we start from the original formulation of the galileon model 
with the original choice for couplings to gravity. Within this framework we find that there 
is still a subset of possible Lagrangians that give selfaccelerating solutions with stable 
spherically symmetric solutions. This is a certain constrained subset of the third order 
galileon which has not been explored before. We develop and explore the background 
cosmological evolution of this model drawing intuition from other even more restricted 
galileon models. The numerical results confirm the presence of selfacceleration, but also 
reveals a possible instability with respect to galileon perturbations. 
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1. Introduction 

Supernova data (l], in conjunction with data from the CMB Q suggests that the Universe 
is expanding at an accelerated rate. Though alternatives to accelerated expansion have been 
proposed in terms of inhomogeneous models (see for instance [||] ) , it is now most commonly 
asked in cosmology not whether, but how or why the Universe's expansion is accelerated. 

The simplest way to get accelerated expansion is by adding a cosmological constant. 
This solution is now part of the standard cosmological model and it fits well with modern 
precision data of the cosmic microwave background, supernova data and baryon acoustic 
oscillations |3]]. 

However simple and elegant the cosmological constant solution may seem, questions 
still remain. If this really is the right explanation for the cosmic acceleration, what is the 
physical origin of the constant? Why has it got the value that it has? A possible physical 
source is vacuum energy, but attempts to predict its value from known microscopic physics 
notoriously yield results more than a hundred orders of magnitude larger than the observed 
values [||. Though there are proposed resolutions to this, appealing both to our lack of 
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knowledge of fundamental physics and to anthropic arguments ||, one might alternatively 
look for a dynamical solution important on cosmological scales. 

The alternatives to the cosmological constant come in roughly two categories. The first 
proposes some sort of generalised dark energy ||. The other possibility is that of some 
kind of modification of gravity on large scales @. In most cases models of dark energy and 
modified gravity are dual in the sense that a conformal transformation can take a theory of 
one type into a theory of the other kind. The main difference lies in what the formulation 
is like in the frame where one assumes the background metric. If the theory in this case is 
minimally coupled to gravity, the theory is considered a dark energy model, whereas if it 
is nonminimally coupled to gravity in this frame it is considered a modified gravity model. 

The latter alternative consists of such alternatives as additions of generalised metric 
functions, /(-R)-gravities || and the Gauss-Bonnet-gravities [||, 10 1, scalar- vector-tensor- 



theories |ll|, |12J, additions due to extra dimensions |l_3[ such as the DGP model |14| , |15] 
and its generalisations, degravitation jl(]] and cascading gravities [17|, the Brans-Dicke 
theories gj| ||, || and many more 



Recently the Galileon model [21] is a modified gravity that has attracted a good deal of 
attention, though it has also been formulated as a pure dark energy candidate. Motivated 



by the DGP and its Vainshtein mechanism, the authors of [21] constructed the Galileon 
model, a model of a scalar field n that has a symmetry of the equations of motion under 
the galilean shift: 

(1.1) 



7T 



7T + C + 



and found all five possible Lagrangians that obey this symmetry in four dimensions. The 
five different terms are denoted by the order of the scalar field it they contain, and there 
is one term unique up to total derivatives of each of the orders one to five, and none of 
higher than fifth order. They considered a theory built from linear combinations of these 
terms and found that it could yield selfaccelerating solutions in a cosmological background. 
By considering spherically symmetric solutions and their stability, bounds were put on the 
coefficients. 

After this promising start considerations of the cosmological evolutions of galileon 



theories or subsets of such theories were considered in [22], which dealt with only the 
behaviour of the third order galileon term, which turned out to give a slight degravitation 
of the cosmological constant, but gave no viable selfaccelerating solutions, and for slightly 
different couplings to gravity than those given in [^] in [^, 24, ^5|] which all explored the 



cosmology of theories of up to fourth and in the case of [£4| fifth order galileon models. 

Several other theoretical aspects of the galileon theory have also been considered, such 
as the null-energy-condition (NEC) violation |2q| , covariantisation [27], the Vainshtein 
mechanism and its effect on spherically symmetric solutions [^] and black hole accretion 

m. 



From theoretical considerations of the symmetries, many generalisations have been 
built. In pO] expressions for all the galileon invariant Lagrangians in arbitrary dimensions 



were derived. Generalisations to arbitrary p-form galileons were obtained in ]31] and the 
multi-galileons, which are multiplets of scalar galileons, contained in this formalism has 
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since been considered in [^] and which mainly focused on the two- or bi-galileon case, 
in 1 34] which focused on theoretical aspects and the evasion of Derrick's theorem [35, 36] 



and [37] which considered instabilities in spherically symmetric solutions. The paper [38] 
showed that the N-galileons are the goldstone bosons of the spontaneous symmetry breaking 
of the N + D dimensional diffeomorphism group in forming a codimension N brane of D 
dimensions. This property of the multigalileons was also mentioned in [32| where it was 
noted that the bi-galileons were the goldstones from cascading gravity from 6 dimensions. 

The fact that the galileon is a natural artifact of a five dimensional theory was further 
reviewed in [39] which found in this a unification of the DBI [40, |4l], 43, 44] and 
galileon theories. The two theories turn out to be different nonlinear realisations of the 
5D diffeomorphisms on a 4D brane. A generalisation of this theory was formulated and 



explored in [45]. 

The connection between massive gravity and higher dimensional theories are well 
known, and in the case of the galileon the connection to massive spin 2 fields was shown 
in pL 



Many other generalisations of the galileon model have also been explored. In [47] 
a generalisation where n — > /(vr) was considered in a cosmological context. Some of 
these generalisations break the galilean invariance. For instance a generalisation to the 
coefficients of the different galileon Lagrangian terms being functions of the galileon field 



was considered for up to third order galileons in | 4q| , |49|, |56f| and |19 |. Up to third order 
galileons have also been generalised to theories with coefficients not only dependent on the 
galileon field itself, but also on its kinetic term X = —g^ v d^ird v iT/2 in ||5l"|] which considered 
perturbation equations in such a generalised scenario, in []52| where cosmological evolution 
was reviewed and in [1531 which considered inflation. 



Though the original formulation of the galileon was proposed as a dark energy alterna- 
tive [|21[ , the model has recently been studied as an alternative to inflation. In particular 
its NEC-violating properties have attracted attention in this context as this can give rise 



to a galileon genesis from flat space [54] yielding a theory with no Big Bang. The NEC- 



violation was also found in the generalised theory of [53] in the inflationary scenario. In [55] 
a covariant version of the galileon is used to drive inflation, motivated by the fact that the 
constant shift symmetry is already only mildly broken in the traditional slow-roll inflation, 
and hence the galileon is a natural and useful place to look for generalised inflation. 



We want to take the model back to its basic formulation in [21], where we feel that a 
simple self-accelerating solution has yet to be properly explored. In this paper we find from 
the theoretical constraints in [21] and [26] that a third order galileon theory with a first 
order term, a so-called tadpole, seems to be allowed, stable and yielding selfaccelerating 
solutions in a certain regime of the parameter space. We explore this simplest of the 
galileon selfaccelerating solution and its background evolution, as we feel that a thorough 
treatment of this theory seems the natural place to go first from the results of the original 
papers [21] and p6|. A treatment of the perturbations for this model will be postponed 



for a future publication. 
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2. The galileon Lagrangian 



Before we start the treatment of our particular galileon model, we would like to recapitulate 
the formulation of the galileon model given in [^TJ . 
Looking at the action: 

S = J <£x^-Ml L ^n\R + ± VT^ + C n + ttT£] (2.1) 

where Mpl is the four-dimensional reduced Planck mass, R is the Einstein frame Ricci 
scalar, T^ u is the stress-energy tensor for the matter fields, \h^ v is the coupling between 
gravity and the matter fields compensating for the matter not feeling the same metric as 
the Einstein frame one and tt is the galileon. The equations of motion of the galileon 
Lagrangian is invariant under the galilean transformation 

tt -> tt + c + (2.2) 

In 4D there are 5 galilean invariant Lagrangians, one for each order of the galileon up 
to 5. Higher order terms involve only total derivative terms. Defining II = d^d u ir, [. . .] to 
be the trace of the operator . . . and • to denote Lorentz-invariant contraction of indices, 
the five galilean invariant Lagrangians are: 

(2.3) 
(2.4) 

(2.5) 

dir - [n 2 ]d7T ■ d-K + 28tt ■ n 2 • 8tt) (2.6) 



+6[n]5vr • n 2 • 8tt + 2[n 3 ]9vr • dir + 3[n 2 ]5vr • II • dir 

-68tt • n 3 • <9vr) (2.7) 

(2.8) 

In general, then, the full galilean field Lagrangian takes the form 

5 

£ 7r = ^c i £ i (2.9) 

i=i 

For our selfaccelerating solution we want a stable deSitter attractor for a vanishing 
matter stress-energy tensor. Considering in this case the effect of the galileon on the 
cosmological evolution yields a configuration for the tt in terms of the wanted asymptotic 
value Hq of the Hubble parameter in this case. Plugging this into the equation of motion 
for the galileon leads to the following condition on the coefficients of the Lagrangian 1 

ci - 2c 2 Hl + 3c 3 H* - 3c 4 Hl + ^c 5 H$ = (2.10) 



A = 


TT 




c 2 = 


— -drr ■ Ott 
2 




c 3 = 


-i[n]5vr-97r 




u = 


~ ([n] 2 9vr • 8tt 


- 2[n]d7r • 


c 5 = 


--([UfdTT-dTT 

5 V 


- 3[n] 2 3vr 



1 For the full derivation of this condition in detail see |2ll 
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Now, we want the internal galilean symmetries and external Lorentz symmetries to be 
unbroken, so the perturbations should obey these symmetries, hence if our equations of 
motion for the galilean terms are given by 



5tt 



(2.11) 



then the perturbations to the equations of motion must be given by 



(2.12) 



where di are the coefficients for the perturbations. To simplify the calculations we assume 
the coefficients d{ to incorporate the zeroth order deSitter solution plus the perturbations. 
We can then find invertible relations between the djS and the as for i = 2 . . . 5. 



C2 




C3 




C 4 




_ C 5_ 





1 3H$ 3H o 



1 





3H$ 3H$ 



2H$ 
1 



ci is then simply given by the equation ( p. 10 ). 
The inverse matrix becomes 





~d 2 ~ 




d 3 




di 




_d 5 _ 



~ d 2 ~ 




ds 




di 




_d 5 _ 





1 - 

1 





o u2 9 ttA 
6N 2 H 



-3H$ 



-3Hl 3H^ 



1 





1 





C2 




C3 




C 4 




. C 5. 



We should now check what constraints must be put on these coefficients for the solution to 
admit a stable deSitter solution with spherical symmetric solutions about compact sources. 
These configurations should be ghost free. We also want the fluctuations to be subluminal. 
Combining all this demands the conditions for the c^s are given by 



d 2 > 
d A > 



d 3 > y -jd 2 d± 
d 5 < 



(2.13) 



We should study the cosmology of such a model where the null-energy condition is vi- 
able, and see whether it can reproduce results for ACDM cosmology in terms of background 
equations and perturbations. 

Actually the constraints above already ensure subluminality and no breaking of the 
NEC. In the coming subsection we check that this is actually true in the case considered 
in the article [26|. 
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2.1 The case of C\ = 

In the article (2(| they consider a deSitter solution with c\ 
of no violation of the NEC conditions means that: 

2 c 2 1 



C4 < 



333 



+ 



f'3 



111F 2 



and check for the demand 



(2.14) 



When ci = 0, then C5 is completely determined from the deSitter criterion ( 2,1C| ). 
This means that we can rewrite our transformation matrices so that they only transform 
between indices 2, 3, 4. The results are 



C2 




C3 




C 4 









ds 




C?4 





-3 


-3i?< 


4 

8 


-5 

4 


ml 




"-3 


3F 2 


4 


-5 



3H„ H; 



i) 



3 tr4 " 
~2 H 




~d 2 


-3HI 






-3 






~2 H 




C2 






C3 


-3 




C 4 



By using the transformation matrices above to turn the condition on cs given in equa- 
tion into a condition on the (is, we realise that the conditions given in equation (|2.13j ) 
are sufficient to prohibit the breaking of the NEC. 



3. A possible Lagrangian for the galilean modified gravity with selfaccel- 
erated solutions 

Since the demands on the coefficients already ensures a deSitter solution, we should look at 
the early evolution and the evolution of perturbations in this theory. Since the assumptions 
also imply stability of the deSitter solution we should only check the stability in the more 
general non-deSitter regime. 

However, constraints from spherical sources compels us to choose d^ = d$ + C4 = C5 = 
0. 2 Hence we are left only with C2 and C3. The constraint in equation ( p. 13 ) becomes 



ci - 2c 2 Hl + 3c 3 H% = ci - 2d 2 Hl - 3d 3 H$ 







(3.1) 



and requires c\ to be larger than zero for the constraints on the ds to hold. Hence we get 
a Lagrangian: 



(3.2) 



where 



C3 > 

c 2 > 3H$c 3 
ci 



2c 2 tf 2 - 3# 4 c 3 



(3.3) 



2 See fell 
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where Hq is the value of the Hubble rate in the deSitter attractor solution. In the later 
numerical computations we set this value to the present day Hubble rate. 

We choose to work in a Jordan frame where matter is minimally coupled to the galileon, 
thus considering the galileon a modification to gravity as was eventually also done in the 
formulation made by 21]. Our final action becomes 

S = J d A x^-g (^k e - 2 *R + Ci tt - | (dir) 2 - |Dvr(9vr) 2 ) + S matter (3.4) 

where Setter is t ne action of other fields like photons, neutrinos, baryons and dark matter. 

Because of the coupling between the galileon and gravity the resulting Einstein equa- 
tions and equation of motion for the galileon field are no longer galilean invariant in a 
non-flat background. However, since the Jordan frame has been chosen a coupling between 
gravity and the galileon must be supplied. This particular choice makes the action reduce 
to the ordinary Einstein-Hilbert action in the limit of no galileon fields, whilst keeping 
the galilean symmetry in flat space. Hence this choice for a coupling to gravity is natu- 



ral. It is also the same non-minimal coupling as chosen in [21] so all their stability and 
selfacceleration conditions hold for this model. 



From this and reviewing the result of 22] we can find the equations of motion. For tt 

it is: 

ci + c 2 Dir + c 3 ({Ott) 2 - (V^V^) 2 - R^V^V„n) = M$ L Re- 2n (3.5) 



Using the result from |£2[ combined with the derivation in [56] the Einstein equations 
become 



(3.6) 



^ (2V M 7rV J/ 7rD7r + ^V a 7rV Q (Stt) 2 - 2V (At 7rV, ) (c>7r)< 



This should enable us to study background cosmology and eventually also perturbation 
theory in this modified version of gravity. 

3.1 Equations of motion in a homogeneous and isotropic Universe 

To simplify calculations and the considerations of the equations of motion derived above we 
start by assuming isotropy and homogeneity. We also assume the Universe to be spatially 
flat. This means that we are working in a flat FRW metric background. This derivation is 



very similar to the one presented in section 4 of [22|. 

In the flat FRW metric background, the equation of motion for tt in equation (ft.5| ) 
now becomes: 



M^ L 6(- + H 2 )e~ 2n = ci - c 2 (# + 3fl7r) + 



a 



3c 3 (I {Hi: 2 ) + 3H 2 n 2 ^) (3.7) 
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where we see that the term proportional to C3 is equivalent to the term found in [Effil . 

The 00 component of the Einstein equation gives the equivalent of the modified Fried- 
mann equation 

3e' 2n M^ L H 2 = p - civr + ^-rr 2 + 6H-kM^ L e- 2w - 3Hc 3 n 3 (3.8) 
The combination — i (00 + yields the Raychaudhuri equation 

M 2 L e" 2 ^ =~(P+ 3P) " I (civr + c 2 tV 2 ) + |* 2 (H-k - n) 

+M$ L e~ 2n (Hit + tt- 2tt 2 ) (3.9) 

which has a slight discrepancy of e~ 2n HirMp^ with respect to the one found in |£2| which 
seems to be due to a difference in calculations of covariant derivatives. 



As in 1 22] we can find an effective energy density p^ and pressure P n for the ir field. 

Pjt = -c l7 r + ^tt 2 + 6HnM^ L e~ 27T - 3Hc 3 ir 3 (3.10) 

P n = Cl vr + ^tV 2 + c 3 t\t\ 2 - 4M£ L Hire- 2n 

-2M| L 7re -27r + 4ikq?. L 7i-V 27r (3.11) 



Again the part of P w resulting from C 3 is not exactly the same as found in |22J. The same 



discrepancy, presumably caused by the same calculations difference, is still present. 
4. Numerical simulation of the background equations 



Before we start the numerical simulations of the background equations ( |3.7D , ( |3.8|) and 
fl3.9|) it is useful to rewrite the derivatives with respect to t in terms of a new variable 



x = In a. Since = H4- we get for the equation of motion of the galileon: 
M^ L 6H [H' + 2H) e~ 2n = c x - c 2 H (flV + Hit" + 3Htt') + 

3c 3 (^(^ M 2 ) +3 ^(vr') 2 ) 

= ci - c 2 H (HV + #V' + 3iZV) + 
3H 3 c 3 tt' (3i2V + 2Ftt" + 3i?V) 

where ' denotes a derivative with respect to x. 

The Friedmann and Raychaudhuri equations come out as 



(4.1) 



and 



3e- 2n M$ L H 2 = p - Cl vr + ^-H 2 (tt') 2 + 6tf V 'M 2 L e~ 27r - 3tf 4 c 3 (tt') 3 (4.2) 

M 2 L e~ 27T H (H' + H)=-~(p + 3P) - i ( Cl vr + c 2 tf 2 (tt') 2 ) + 

|tf 3 (tt') 2 (#7^ - tfY - Fvr") 

+flMp L e _27r (#V + flV + Hit" - 2H (tt') 

(4.3) 
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We now take advantage of the nice property of the Jordan frame that matter is not 
coupled to the galileon field. This means that the equations that determine the development 
of the matter or radiation from their equations of state are not altered, so p m = f2 m a~ 3 = 
Q m exp (— 3x) and p 7 = Q 7 a -4 = f2 7 exp (— 4x), p = p m + p 1 and p + 3P = p m + 2p 7 . 

We use the above equations to make a set of three coupled differential equations to 
solve numerically in H, tt and vr'. First we use the equation: 



thereby transforming the set of second order equations ( ^4 . l| - |4 . 3| ) in vr into a set of first order 
equations in vr, tt' and H. Then we need an equation for vr" and one for H' . However, we 
need at least one of these not to depend on the other. 

We enter the result from the first Friedmann equation into the Raychaudhuri equation, 
to simplify slightly. The result is: 

Mhe- 2 *H (h> + *H)=-lp- \c x tt - \c 2 H 2 (tt') 2 - f (vr') 2 (W + Htt") 

+HM^ L e~ 2n (2Htt' + H'tt' + Htt" - 2H (tt') 2 ' 



(4.5) 

Then using ( f4.5|) to solve for the second derivative of the galileon field we get: 

vr" = 7 + 5H' (4.6) 



where 



Mp L H 2 e~ 27T (1-2^ + 2 (ir') 2 ) + \P + \c x tt + \c 2 H 2 (vr') 

T = 

H 2 (M| L e-2^ _ ^H 2 (tt')'' 

. M 2 L e~ 2 * (1 - vr') + f H 2 (vr') 3 



H [M 2 ^-^ - f H 2 (vr') 
Plugging this into equation (|4.1|) we get that 



(4.7) 



„,_<*- c 2 H 2 (7 + 3vr') + 3# 4 c 3 vr' (2 7 + 3vr') - 12H 2 M 2 L e- 2n 
6HM 2 L e~ 27r + c 2 H (vr' + H6) - 3c 3 ff 3 vr' (3vr' + 2H5) 



(4.8) 



For reference it is nice to have the formulas for the galileon density and pressure given 
in terms of the ' derivatives: 



_ Cl7r + C A H 2 (^'f + 6H 2 7T>M 2 L e- 2 * 



3# 4 c 3 (vr') : 



P„ = Ci tt + ^H 2 (tt'Y + c 3 H 3 (vr')" (H'tt' + Htt") - AM 2 L H 2 Tr'e 
-2M£ L He~ 2n (H'tt' + Htt") + 4M| L # 2 (vr') 2 e~ 27T 



2 / -2tt 



(4.9) 



(4.10) 
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5. Results in different particular regimes 



Before moving on to our full model it is interesting to review and study the behaviour in 
different regimes of special values for the coefficients. Some of these may be interesting 
for their simplicity. Others may be comparable to theories considered elsewhere. Together 
they provide useful intuition for our interpretation of the evolution of the full model. 

5.1 All coefficients zero 

If we set ci = C2 = C3 = 0, the galileon field will still be present and dynamical, and will 
have a special effect on the evolution of the Universe. We can see this directly from writing 
out the equations in this case: 

The galileon equation of motion gives: 



- + H z = 

a 



(5.1) 



Adding the Friedmann equation and Einstein equations together and writing p tot = p + p w , 
Ptot = P + Pn, we get: 



- + H 2 = -- tot + 3P to t) + 7: Ptot, 
a 6 3 

which leads to a demand for an effective total equation of state at all times of 

Ptot 1 



W e fi 



Ptot 3 



(5.2) 



(5.3) 



This means that the Universe appears radiation dominated at all times. Numerical results 
show this, as can be seen from the evolution of H shown in figure |]. 




Figure 1: Plot of the evolution of the Hubble parameter H in the galileon model with all coefficients 
ci = c 2 = c 3 = 0. Note that H evolves as if always in a radiation dominated Universe. 
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5.2 Tadpole only, C2 = C3 = 0, non zero c\. 

The case of c\ being the only non-zero constant greatly resembles the case considered 



above in section |5.1| with no coefficients. However we now have a linear potential making 
the theory start off in such a way that the evolution starts not like a radiation dominated 
one, but one defined by the constant c\ if the galileon starts out very small. We see this 
from the equation of state of the galileon in this case: 

Sl] +H 2 = ^-^L- (5.4) 

As the galileon grows the term on the right will grow leading eventually to accelerated 
expansion. If the galileon field on the other hand has negative values, the evolution will 
move towards the radiation like one as the absolute value of the field becomes larger. 
Adding together the Friedmann and Raychaudhuri equations we get that 

c\ = ptot (1 - 3w eS ) (5.5) 

We note that as c\ 7^ the theory can only approach a radiation dominated like epoch 
asymptotically as p t ot gets very large. 

In the tadpole only case the equations for the numerical simulations become: 

7 = - - 2vr' + 2 (vr') + 



2 v ' 2H 2 M^ L e-^ 

1 - 7r' , , 



and H' is simply given by equation (5.4). 

For very small values of c\ we of course retain results close to those in the case discussed 
in section 5.1. For larger values of c\ the results are shown in figure |2] for the evolution of 



H and in figure ||] for the evolution of the equation of state parameter. 

From the two plots we see that in the beginning the Universe behaves as if radiation 
dominated, or very nearly so. This is because early on pt ot is very large so that we approach 
the asymptotic value for the equation of state parameter w e s = 1/3. As the evolution goes 
on it deviates from this and starts showing a smaller equation of state parameter and a 
smaller deceleration. This shift happens earlier, giving smaller values for the equation of 
state parameter and slower deceleration the larger the value of c\. For sufficiently large 
values of c\ a period of acceleration can even be reached. However, after experiencing a 
minimum value of the equation of state parameter or deceleration, the equation of state 
parameter starts growing again eventually asymptotically towards an evolution with < 
w e g < 1/3. We see that by making c\ large enough we can get a value for H today equal 
to the one we find from current observations 3 . However, in these cases the accelerated 
expansion has already taken place and passed, the Universe now being in a decelerated 
state. We also see that most of the evolution history will change considerably in this case. 

Considering specifically the two cases of largest c%, c\ = 10 6 Mp L /?Q and c\ = 10 8 Mp L f^Q 
since these have the most interesting evolution, we plot the values of the different densities 
and fractional densities in figures ||, ^ |6| and [/]. 
3 See for instance 10]. 
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Figure 2: Plot of the evolution of the Hubble parameter H in the galileon model with coefficients 
C2 = C3 = and nozero c x . The evolution of H is shown for five different values of c x : c x = 
0.01M| L # 2 , ci = lOOM^Hg, ci = 10 4 M| L i^, c x = 10 6 M^ L H§ and c a = 10 8 M| L JI|. The 
evolution of iJ for the usual ACDM model has also been plotted alongside the others. 
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Figure 3: Plot of the evolution of the total equation of state parameter w e ff in the galileon 
model with coefficients C2 = C3 = 0, c x nonzero. The evolution is shown for six different values 
of c x : ci = O.OlM^Hfi, ci = Mp L i?o) c x = 100M| L i?£, c x = 10 4 M| L .ff^ c x = lO e M* L H$ and 
ci 



5.3 c x 



C3 



0, non zero C2 - the Brans-Dicke type theory 



The case of only C2 different from zero is in fact a type of Brans-Dicke theory [57]. We can 
see this from starting with the generic Brans-Dicke action: 



Sbd 



(j)R 



(5.7) 
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Figure 4: Plot of the evolution of the fractional densities Qi in the galileon model with coefficients 
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Figure 5: Plot of the evolution of the fractional densities fli in the galileon model with coefficients 

-2 tt2 



c 2 = c 3 = 0, ci = f0 8 Af| LJ ff^ 



To get this into a form similar to our own we realise that we must set: 



(5.8) 



and inserting this we see that for the second term to correspond to our C2 term we must 
have: 

(5.9) 



I.e. a theory with only the C2 term present is equivalent to a Brans-Dicke theory with the 
Brans-Dicke parameter lo(4>) oc 4- 



- 13 - 




le+05 10000 



Figure 6: Plot of the evolution of the densities pi in the galileon model with coefficients ci = C3 = 0, 
ci = 10 6 Afp L i?Q. Since the galileon density starts out negative its absolute value \p v \ is shown in 
the plot. 
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Figure 7: Plot of the evolution of the densities /?,: in the galileon model with coefficients ci = C3 = 0, 
ci = f0 8 Afp L i?Q. Since the galileon density starts out negative its absolute value |p w | is shown in 
the plot. 



The background equations we need to simulate now have: 



1 - it' 



2 4H 2 M£ L e~ 2 ' 



H < 5 ' 10) 
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Plugging this into equation (|4.1| ) we get that 

tti = c 2 g 2 ( 7 + 37rQ + 12ff 2 Mg L e- 2 - 

6HM^ L e-^ + c 2 H(tt' + H5) { ' ; 

The results for the evolution of the Hubble parameter H and the effective equation of 
state parameter for different values of C2 is shown in figures |8| and ^ respectively. 




Figure 8: Plot of the evolution of the Hubble parameter H in the galileon model with coefficients 
ci = C3 = and nozero c 2 . The evolution of H is shown for five different values of c 2 : c 2 = 0.01Mp L , 
c 2 = 1M| L , c 2 = 10M^ L , c 2 = 50Ai| L and c 2 = 10 8 M^ L . The evolution of H for the usual ACDM 
model has also been plotted alongside the others. 




Figure 9: Plot of the evolution of the total equation of state parameter u> ff in the galileon model 
with coefficients c\ — C3 = and nozero c-i. The evolution is shown for five different values of ci'- 
c 2 = 0.01M! l , c 2 = 1M| L , c 2 = 10M| L , c 2 = 50M| L and c 2 = 10 8 A/| L . 
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We see that C2 in fact acts as a demarcation scale between the non-galileon background 
theory and the theory with a galileon with all the coefficients equal to zero. The larger 
the value of C2, the later the galileon will have any impact on the evolution. However, as 
we raise the initial values of the galileon, the larger is the C2 we need to compensate for 
this. Hence quite large initial values for the galileon will make the galileon dominate and 
give a Universe at early times with much smaller values of the Hubble parameter before 
transitioning into the radiation like phase. For instance in the case of C2 = 50Mp L the 
initial values of ir and ir' were set to 0.5e — 10 and 0.5e — 50 respectively in the case shown 
in figures |8] and ||. Changing the values to 0.5e — 10 and 0.5e — 10 respectively had very little 
effect, however a change to 0.5e — 5 and 0.5e — 5 for the initial values respectively changed 
the evolution of the background quite considerably, resulting in a much lower value for the 
Hubble constant, reaching a value comparable to that in the plotted case for c 2 = 0.01Mf L 
only very recently (redshift of only a small fraction), and from there transitioning into 
radiation like behaviour. 

Since the theory only shows transitions between a pure radiation like model and the 
theory without galileons, this theory can not display selfacceleration. Hence for it to give 
accelerated expansion a cosmological constant or some other form of dark energy would 
need to be added. 

That the impact of the galileon decreases for increasing values of C2 can also be seen 
by comparing plots of the fractional densities and densities for two different values of c^- 
This is shown for C2 = 10Afp L and 10 8 M| L in figures FRl O, O and O respectively. 
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Figure 10: Plot of the evolution of the fractional densities fli in the galileon model with coefficients 
d = c 3 = 0, c 2 = 10A^ L . 



5.4 Non-zero value of C3 and added cosmological constant 



The theory considered in [22] corresponds to a model of non-zero C3, c\ = C2 = with the 
addition of a cosmological constant. That is, a theory of a third order only galileon in a 
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Figure 11: Plot of the evolution of the fractional densities fli in the galileon model with coefficients 
d = c 3 = 0, c 2 = 10 8 M2 L . 
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Figure 12: Plot of the evolution of the densities pi in the galileon model with coefficients c\ = 
C3 = 0, C2 = fOMp L . Since the galileon density starts out negative its absolute value \p v \ is shown 
in the plot. 



ACDM background. 

The equations in this case simplify significantly. The equations for the coefficients 7 
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Figure 13: Plot of the evolution of the densities pi in the galileon model with coefficients c\ = 
C3 = 0, C2 = 10 8 Afp L . Since the galileon density starts out negative its absolute value |p w | is shown 
in the plot. 



and 5 given in equation (4.7) become: 
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1 - TT') + M 2 
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7T 



if [ M2 L e- 2 ^ 



f {TT'f 



and the equation for H' becomes: 

H 3 c 3 tt' (2 7 + 3vr') 



4i?M2 L e- 27r 



2M2 e- 2 ^ 



(5.12) 



(5.13) 



M>L . - c 3 HV (Sir' + 2HS) 

We can obtain a model equal to the one considered in p2| by considering this regime 
with ci = C2 = 0, C3 7^ in a background with a cosmological constant. However, since we 
had a slight discrepancy in the equations as compared to their results, there might be some 
differences. We must also pay attention to the units, as we are working with dimensionless 
7r fields, where as they have dimensionful tt fields. 

To review the results and compare to those of [^] , we have plotted the evolution of the 
Hubble parameter H as compared to that of ACDM, the evolution of the energy densities 



Pi and the evolution of the fractional densities f2j in figures 14, 15 and in figure 16 for a 
particular value of C3, C3 = 100Mp L / Hq. As in [22| we have defined 

SU = j£L • (5.14) 



3H 2 M 2 L 



To make the results even more comparable to those in [22] we make the value such 
that Q m today is close to the value measured today. We also extend the simulation to 
z = -0.8. 
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Figure 14: Plot of the evolution of the Hubble parameter H in the galileon model with coefficients 
ci = C2 = 0. C3 = 100Mp L / Hq. There is also a cosmological constant present in this model set so 
as to ensure £l m today of about the same size as measured. This corresponds quite closely to the 
case considered in |2^[ . The plot also shows the evolution of the Hubble parameter in the ACDM 
case. 
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Figure 15: Plot of the evolution of the energy densities of the different components in the galileon 
model with coefficients c\ = C2 = 0, C3 = 100A/p L /i?g and a cosmological constant. 



The results resemble closely the ones found in [22]. Similarly to [22] we also found that 
a correct evolution history was first obtained for relatively large values of C3, the models 
resembling ACDM more the larger the value of C3. Plots showing the evolution of the 
Hubble parameter and the evolution of the effective equation of state parameter in some 
different cases are given in figures 17 and [E3. These plots show only small differences. 



In figures 19 and EG the difference in evolution histories is shown for even smaller values 
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Figure 16: Plot of the evolution of the fractional energy densities of the different components in 
the galileon model with coefficients c\ — C2 = 0, C3 = 100A/p L / Hq. Here the added cosmological 
constant has been set so as to ensure TO today close to its measured value, as in [^2). 

of C3. In these plots the cosmological parameter has, however, only been set to its ACDM 
value, as normalising it in the very small C3 cases is quite hard due to the large impact of 
the galileon field. 

The tendency is that for small enough values of C3 the galileon field dominates the 
evolution too early for the Universe ever to experience a period of accelerated expansion. 
This should not come as a surprise after considering the case of all the coefficients set to 
zero in the previous subsection. In this case we have seen that the Universe will behave 
as if radiation dominated at all times no matter what the matter content of the Universe 
may be. The fact that the deviation from ACDM is larger with smaller C3 is hence quite 
understandable, since C3 is a demarcation scale, under which the physics still behaves as 
in a non-galileon Universe. The role of C3 is hence quite similar to that of C2 as described 
in section [5J^ and we would expect the model with C2 only with a cosmological constant 
added to have results somewhat comparable to the results found here. 

To check whether the behaviour is altered due to our slightly different results for the 
Raychaudhuri equation ( |3.9| ), we have performed simulations using their equation set as 
well as our own. They yield only slightly different results, i.e. the general behaviour is not 
changed and only very subtle changes in the evolution can be observed by reviewing the 
data very carefully. 

We have also considered the effect of altering the initial values of the galileon field. 
Actually, as long as the initial galileon field is small and positive 7r s t a rt < 0.1, and the 
initial galileon field prime derivative is quite small and negative > vTg tart > — le — 8, the 
evolution histories found for each value of C3 is quite stable, that is with the initial values 
within these bounds, their exact values have little impact on the ensuing evolution history 
of the Universe. Relatively small initial values of the galileon field and its derivative seem 
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Figure 17: Plot of the evolution of the Hubble parameter H in the galileon model with coefficients 
c\ = C2 = 0. C3 is non zero and there is also a cosmological constant present in this model. Three 
cases are shown C3 = 50Mp L /H§, C3 = 100Mp L / Hq and c 3 = 500Mp L /Hq. For comparison, the 
evolution of the Hubble parameter in the usual ACDM model is also shown in the plot. 




Figure 18: Plot of the evolution of the total effective equation of state parameter wtot in the 
galileon model with coefficients c± = ci = 0. C3 is non zero and there is also a cosmological 
constant present in this model. Three cases are shown C3 = 50Mp L / Hq, C3 = 100Mp L / Hq and 
c 3 = 500A4JHI 



like a reasonable starting point, so assuming this model to have the described behaviour 
seems fairly acceptable. 



6. Numerical results for the full model background equations 

Moving on to the full model, with ci, C2 and C3 obeying the constraints of equation (|3.3|), 
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Figure 19: Plot of the evolution of the Hubble parameter H in the galileon model with coefficients 
ci = C2 = 0. C3 is non zero and there is also a cosmological constant present in this model. 
This corresponds quite closely to the case considered in f2l| , however, they have set the value 
of the cosmological constant so as to ensure a flat Universe today. Here I have simply used the 
cosmological constant value of the usual ACDM model. Three cases are shown: c 3 = 0.01Mp L /-ffg , 
C3 = Mp L /HQ and C3 = 100Mp L /HQ. For comparison, the evolution of the Hubble parameter in 
the usual ACDM model is also shown in the plot. 



which is supposed to show selfacceleration, we obtain an evolution for the Hubble parameter 



H as shown in figure 21. We see that the evolution of the Hubble parameter in this case 
is quite comparable to that of the ACDM model. 

However, the evolution in this theory is not exactly as that of a ACDM model. This 
can be seen for instance from the evolution of the densities and fractional densities of the 



different components as shown in figures 22 and 23, or in the evolution of the equation of 
state parameter shown in figure 24. 4 This interesting behaviour may lead to possible ways 



of discerning this model from the ACDM model or at least to possibilities for constraining 
the values of C2 and C3 from distance measurements, such as supernova data or large scale 
structure surveys showing the evolution of the equation of state parameters over time. 

6.1 Robustness of results 

Not surprisingly, the model does not stably yield near ACDM results under all initial 
values of the galileon field tt and its derivative. If these have values that are too large in 
the beginning, the early Universe evolution will be dominated by the galileon and hence 
not give results compatible with current knowledge from baryon nucleosynthesis etc. 

However, the evolution of the background equations also exhibit a dependence on the 
exact values of the coefficients C2 and C3. If C2 is too small, and this seems to not only 
depend on the relative value with respect to C3, the selfacceleration starts too late in the 
evolution of the Universe. If C2 is too large the selfacceleration starts too early. The fine 



The fractional densities are still defined as in equation (5.14) 



- 22 - 







c3 = 0.01 






■ ■ ■ c3 = I 

c3 s 100 








\ 

\ 

\ 

\ ~ 
\ 










- 


10000 


100 

z+I 





Figure 20: Plot of the evolution of the total effective equation of state parameter wtot in the 
galileon model with coefficients c± = C2 = 0. C3 is non zero and there is also a cosmological constant 
present in this model. This corresponds quite closely to the case considered in [p2| . however, they 
have set the value of the cosmological constant so as to ensure a flat Universe today. Here I have 
simply used the cosmological constant value of the usual ACDM model. Three cases are shown: 



c 3 = QMM$JH§, c 3 = MlJHl and c 3 
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Figure 21: Plot of the evolution of the Hubble parameter H in the case where C3 = 3Mp L / H§ , 
c-2 = 18Afp L compared to the evolution of the same parameter in the ACDM case. We see that the 
two evolution histories are quite comparable. 



tuning needed for C2 increases as C3 gets much smaller or larger than a few Mp L /ffg. We 
can summarise by saying that for the cosmological evolution to be close to our own, the 
values of the parameters need to be given by order a few times appropriate factors of Hq 
and Mpl. This makes it tempting to say that although the parameter space is constrained, 
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the allowed values constitute a natural parameter regime. 

After considering the galileon model with each term separately in sections 5.2 to [5.4| this 
is quite understandable. Since C2 and C3 are demarcation scales, keeping the theory away 
from the pure radiation dominated like Universe, at least one of these must be not too small 
to yield sensible results. Since our stability conditions dictate that C2 > 3HqC3, we get that 
C2 must be large enough. Since c\ is the only term that can in fact yield selfacceleration, this 
term must be present, and to get the right value of the Hubble parameter today, and the 
right evolution, c\ must be set by the demarcation scales. Having the demarcation scales 
both too large also causes trouble, since c\ then gets very large and starts the accelerated 
expansion too early. 
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Figure 24: Plot of the evolution of the total equation of state parameter w e s in the case where 
C3 = 3Afp L /i?Q, C2 = 18A/p L . We see that after the radiation and matter dominated epochs, we 
get a slight increase in the effective equation of state parameter before it dips down to negative 
values. After this it starts increasing again. We can interpret the plot as follows. At early times the 
evolution follows that of the standard ACDM model, starting with radiation domination and then 
going into matter domination. As the galileon starts to dominate the evolution it has a negative 
density and is dominated by the coupling to gravity hence driving the theory towards a radiation 
dominated like epoch like in the all coefficients zero case. After only a short while of this however, 
the Ci term takes over leading to a period of accelerated expansion. As time continues to go by 
we cross the demarcation scale given by the C2 and C3 entering an epoch where selfacceleration will 
stop and eventually turn to deceleration again. 

For the initial values of the galileon and its derivative we find that they must be tuned 
to be quite close to zero. In the case of C3 = 3Mp L / Hq, C2 = 18Mp L , the evolution remains 
stable as long as the initial value for tt' < le — 20 with the initial value of ir < le — 9. As 
we change only the initial value of tt to around le — 8, the evolution shortens its matter 
dominated epoch and gets a more abrupt change to the selfaccelerated phase. Outside this 
regime the evolution becomes dominated by the galileon from a very early stage, and it 
becomes totally unstable. This is similar to the results found for perturbations in p0[ . 
where the initial values for the galileon had to be tuned to zero in the beginning in order 
for the galileon fluctuations not to take over and spoil the homogeneity and isotropy of the 
background. 

7. Discussion 

In this paper we have explored a certain version of the galileon model formulated in . 
Using their derivation of stability and selfaccelerated expansion, we have found a model 
that both admits selfacceleration and gives stable spherically symmetric solutions, that is 
a third order galileon with certain constraints on the three parameters. This model has 
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not been thoroughly studied in the literature, and we feel that it is important to provide 
such a study of one of the simplest proposals for a selfaccelerating galileon. 

We have explored the background cosmology of this model, using numerical simulations 
of each term separately to build intuition. In this way we found that it is the tadpole term 
that leads to the selfacceleration of the solution. The second and third order terms defines a 
sort of demarcation scale above which the selfacceleration slows down, creating an evolution 
history not totally equal to that of ACDM. We also found that the galileon density becomes 
negative before it starts dominating. As it passes through zero and on to positive values 
the evolution of the equation of state goes through a phase where it grows slightly. This is 
a feature of the chosen coupling to gravity, which we have seen yields radiation dominated 
like solutions if present with no other galileon dynamics. 

We also found that the parameter space for a solution comparable to ACDM is limited, 
but not totally constrained. The bounds on the parameters to give a ACDM like evolution 
history are, however, consistent with solar system and galaxy cluster gravitational bounds 
found in [21] and they lie in a somewhat natural range. We presume that further studies on 
the perturbations of this model will set even more stringent bounds on these parameters. 

We also found that the initial values of the galileon field and its derivative must be 
tuned to values quite close to zero. With larger values of these initial field, the model 
becomes highly unstable as the galileon begins to dominate quite early, leading to cosmo- 
logical histories very different from the one favoured by present day observations. This is 
similar to results found in |5(]] where the galileon perturbations had to be tuned to exactly 
zero in the beginning to avoid a total destruction of the homogeneity and isotropy of the 
background space. Further studies on the perturbations of this model will be interesting 
in finding out whether the same kind of anisotropic stress will be found here. 

The growth of structure in the galileon Universe will also be of interest as it might 
provide nice signatures in both baryonic acoustic oscillations and CMB, and these can only 
be found through studies of the perturbation equations of this model, which we postpone 
for a future publication. 
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